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The dynamics of a spin- 1/2 particle coupled to a nuclear spin bath through an isotropic Heisenberg 
interaction is studied, as a model for the spin decoherence in quantum dots. The time-dependent 
polarization of the central spin is calculated as a function of the bath-spin distribution and the 
polarizations of the initial bath state. For short times, the polarization of the central spin shows 
a gaussian decay, and at later times it revives displaying nonmonotonic time dependence. The 
decoherence time scale depends on moments of the bath-spin distribuition, and also on the polar- 
ization strengths in various bath-spin channels. The bath polarizations have a tendency to increase 
the decoherence time scale. The effective dynamics of the central spin polarization is shown to be 
described by a master equation with non-markovian features. 



I. INTRODUCTION 

Decoherence is ubiquitous in quantum systems, either 
because of the interaction with the environment Q], or 
because of interventions from the measuring apparatus 
■ An understanding of the phenomenon of decoherence 
is, therefore, essential - from the view point of the foun- 
dational issues as well as the dynamics of open quantum 
systems. In this context, ./V-level systems are particularly 
interesting; the state is defined in a finite dimensional 
Hilbert space ( the associated phase space is compact), 
which makes the theoretical analysis simpler. More im- 
portantly, there is an abundance of N-level systems as 
encountered in NMR systems Q, NQR systems, and po- 
larized photons. In addition, one also has the so called 
pseudo ./V-level systems - occurring in atomic, semicon- 
ductor and quantum dot environments. Although the 
states of the latter systems are defined in infinite dimen- 
sional Hilbert spaces, energetics (at sufficiently low tem- 
peratures) effectively restrict the states to lie in a finite 
dimensional subspace. It is no surprise that the dynam- 
ics of iV-level systems (spin systems in short) has been 
extensively studied. 

Decoherence in spin systems has acquired a new im- 
portance in view of the rapid developments which are 
taking place in quantum information theory. Control- 
lable quantum gates are central to quantum computa- 
tion, and are built of qubits which are physically realized 
through spin-half (or pseudo spin-half) systems. Propos- 
als for qubits include the spin of electrons 0, 0| , nuclear 
spins |(|0and squids Photons are also poten- 

tial candidates since highly entangled states have been 
prepared experimentally by employing the polarization 
degree of freedom Hll Recently, fullerene-based single- 
electron transistors [l^l have also been proposed where 
the electron spin acts as the qubit. Finally, we men- 
tion yet another important example involving the spin 
of the electrons, viz, quantum dot quantum computers 
(QDQC)[b|. An additional feature of QDQC is that the 
number of qubits (electrons) can be controlled precisely: 
starting from zero, electrons can be added one by one. 
The preparation of quantum gates and their manipula- 
tion should be simpler in QDQC. 



Minimising the decoherence involving the electron spin 
in various environments is, therefore, of great impor- 
tance. We study in this paper the dynamics of deco- 
herence with a special focus on QDQC. However, our 
method possesses wider applicability to instances such 
as decoherence in fullerenes, and to spin systems where 
a given reference spin - the central spin - decoheres 
due to its interaction with the bath constituted by the 
other spins. In QDQC systems, the environment is con- 
stituted by GaAs, (Ga, Al)As, and InAs. The nuclei 
69 Ga(|, + 2.016), 71 Ga(§, + 2.562), 27 AZ(§, +3.6414) 
and 115 /n(|, +5.534) have non-vanishing spins and mag- 
netic moments, as shown in the parentheses. The mag- 
netic moments are written in units of nucleon Bohr mag- 
neton. Since the qubit itself carries a large magnetic 
moment, the dominant contribution to the decoherence 
comes from the interaction between the spins. Indeed, 
decoherence times have been measured experimentally 
for qubits both in quantum dot ensembles and single 
quantum dots. These studies find a long decoherence 
time (~ 100ns) in quantum dot ensembles 1141 . Simi- 
larly, studies on a single quantum dot 0, 0, Il6j| indicate 
that the spin flip rate is low, implying a long decoherence 
time again. 

Theoretically, there have been a large number of inves- 
tigations, summarized in a recent paper by Schliemann 
et al 01- A common feature that underlies these inves- 
tigations is the assumption of a mean magnetic field B 
that is produced by the nuclei which constitute the envi- 
ronment. This assertion is based on the observation that 
there are 10 — 10 6 nuclei, all with non-vanishing mag- 
netic moments, that interact with the qubit |l£l llflj. In 
this semiclassical approach, it is the time dependence in 
B that causes the decoherence in the spin state. Reliable 
estimates for the decoherence rates have been obtained 
both analytically and numerically. 

This paper revisits the problem of decoherence with 
a rather different perspective. We show that given the 
phenomenology of QDQC, it is possible to map the hy- 
pcrfine interaction between the qubit and the nuclei 
to an isotropic Heisenberg model with a global inter- 
action between the qubit spin and the total nuclear 
spin. More precisely, the qubit interacts with a spin 
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bath, whose initial state is described by a density ma- 
trix Pb — •Via PIb ■ where the coefficients A/ B are the 
weights for the bath to be in a state with a total spin Is ■ 
Within each spin sector of the bath, the hyperfine interac- 
tion between the bath spin and the qubit spin gets recast 
into an isotropic Heisenberg interaction. With this sim- 
plification, the dynamics becomes exactly solvable: the 
state of the qubit+bath can be determined in a closed 
form explicitly. A partial trace yields the reduced den- 
sity matrix of the qubit. The decoherence rate is easily 
inferred thereafter. 

In the next section we obtain an effective Hamiltonian 
from the basic hyperfine interaction (for QDQC), and 
elaborate on the model in some detail. We then illus- 
trate the dynamics in the simplest of the cases, viz, of 
two qubits in section III. This example is of practical 
importance, exhibiting the so called interaction induced 
entanglement. After the warm up, the general problem of 
a qubit interacting with a spin bath (as described above) 
will be addressed. The expression for the reduced density 
matrix and the dynamics of decoherence will be investi- 
gated at length. Section IV is devoted to examples and 
special cases. Section V shows how the results can be 
extended to more general local Heisenberg interactions 
perturbatively, at small times. In Section VI, we shall 
set up a master equation satisfied by the reduced density 
matrix. The non-markovian nature of the evolution, the 
time-dependent decoherence rates, and the unitary part 
of the evolution will be clearly identified. 



II. THE MODEL 

A dominant interaction in QDQC between the electron 
spin and the nuclei is given by the hyperfine coupling of 
the electron spin to the nuclear spins: 

H hf = Y / K i S-I i (1) 

i 

where Ki is the interaction strength between the qubit 
spin S and the i'th nuclear spin Ii. The coupling con- 
stant Ki depends on the basic coupling strength and the 
wave function of the electron, Ki — K \ ^ (f) |p. =ff . , where 
the electron density is evaluated at the location of the 
nucleus. Merkulov et al 20] have pointed out that the 
spin orbit interaction, which is otherwise dominant in 
spin decoherence, is suppressed here because the qubit is 
strongly localized. For the same reason the contribution 
of the phonon mediated interactions is also suppressed 
at low temperatures [2*oj . Thus the hyperfine interac- 
tion between the electron and the nuclei emerges as the 
dominant mechanism for spin relaxation. 

In addition to the hyperfine interaction, the electron 
spin may also be subjected to an external magnetic field 
which is responsible for its polarization. Alternatively, 
a spin polarized electron can be directly injected into a 
quantum dot, as Cortez et al|2l| have done. Remarkably, 



they have also succeeded in 'writing' and then subse- 
quently 'reading' a quantum dot. We shall be interested 
in this "zero field" scenario. In either case, the nuclear 
coupling to the field is relatively negligible because of the 
smaller values of the nuclear magnetic moments. Finally, 
mention should be made of the internuclear dipolar cou- 
pling. This coupling (strength is of order 10~ 12 eV) is a 
small perturbatio to the other interaction terms. But it 
has an important role in determining the initial spin state 
of the nuclei, and its effect on the subsequent dynamics 
is negligible [l3|. 

As mentioned earlier, each electron is in an environ- 
ment of 10 4 — 10 6 nuclei. It has been assumed, there- 
fore, that the cumulative effect of the nuclear spins is to 
produce an effective magnetic field B - the Overhauser 
field - at the site of the qubit. Likewise, the action of 
the electron spin on the nuclei makes B time dependent, 
which in turn causes spin flip transitions in the qubit 
(l7L llsl |20| . For an explicit determination of the tran- 
sition amplitudes, computational schemes in this semi- 
classical approach involve either taking a gaussian distri- 
bution for B Il8| or considering special configurations for 
the initial state |22ll23| . The contribution of intra-nuclear 
interaction to the time-dependent magnetic field, and the 
resulting spectral diffusion of the qubit polarization were 
considered by de Sousa and Das Sarma0, and Deng 
and Hu 2[| . Let us summarize the semiclassical results 
in brief. Of the two decoherence times, it is is found that 
the phase decoherence time is larger; standard spin flip 
transition calculations (using the Fermi golden rule) yield 
decay profiles which are nonexponential [23| . More inter- 
estingly, the same authors claim that for a large Zeeman 
field, the polarization of the electron decays ~ t^^ . This 
surprising result is obtained for a special choice for the 
initial state of the qubit-nuclear system. Numerical sim- 
ulations by Schliemann et al 19] for a variety of nuclear 
spin states have yielded several interesting results: (i) the 
polarization of the electron decreases in magnitude, with 
accompanying oscillations whose period is estimated to 
be Att/K, with only a weak dependence on the initial 
state; (ii) if the nuclei are initially only randomly corre- 
lated, the spin dynamics of the electron is significantly 
different in that the decay is much faster. A master equa- 
tion for the qubit density matrix was set up, and the 
dynamics was investigated perturbatively by Coish and 
Loss [2^, where the effect of the electronic wave function 
on the qubit decoherence was also considered. 

Interesting that these results are, the semiclassical ap- 
proximation needs a validation from a more rigorous 
analysis. Indeed, the hyperfine Hamiltonian is itself not 
easily amenable to a full quantum treatment. We never- 
theless argue that the phenomenology of QDQC allows 
us to map the problem to the more tractable isotropic 
Heisenberg interaction. To that end, consider a circular 
geo metry for the quantum dot containing the electron 
(l9j . Taking the confinement to be either parabolic or 
coulombic, the electron wave function ip(r) is either a 
gaussian or an exponential with the distance. The ef- 



3 



fective coupling to the nuclear spins which are farther is 
accordingly suppressed, with the maximum contribution 
coming from the nearest nuclei, all of which are roughly 
equidistant from the electron. Indeed, in GaAs there 
are about 45 nuclei in a volume of 1 nm 3 , and the elec- 
tron wave function is roughly uniform over a distance 
of 2-3 nm (wh ile the size of the quantum dot is about 
20 nm) fvA Wl\ . This translates into about a few hun- 
dred nuclei that interact with the qubit with the same 
coupling strength Ki in Eq^ and thus the Hamiltonian 
assumes a simple form 

H eff = KS-T B , (2) 

where Ib = ^Zh is the total spin of the (nearest) nuclei. 
If we were to consider the nuclei in the next nearest cir- 
cle (next to nearest neighbors in these two dimensional 
samples) , the Hamiltonian acquires additional terms with 
another effective coupling constant K' whose value is ex- 
ponentially suppressed relative to K, and may be treated 
as a (small) perturbation. We shall ignore all such higher 
order contributions here. The dynamics has been ef- 
fectively mapped, to an excellent approximation, to an 
isotropic Heisenberg interaction between the qubit spin 
and the total nuclear bath spin. 

The above Hamiltonian can have a more general role: 
consider an assembly of N identical spin-half particles 
interacting through a Heisenberg interaction. We can 
write the Hamiltonian of the spin system as 

Hspin — ^ ^ KijSi ■ 5j , (3) 

where Kij is the interaction strength of a pair of spins Si 
and Sj. For any chosen spin, denote it by S, the other 
spins would constitute a bath. We can split the above 
into terms involving the interaction of the spin and the 
bath, and an intra-bath term, as 

H spin = K X S ■ I Bl + K 2 S ■ I B2 + H bath , (4) 

where we displayed the interaction strength K\ (K2) of 
the central spin with the total spin of the first (second) 
neighbours Ib 1 {Ib 2 ) ■ And similarly there are interactions 
with further neighbours. The last term is the intra-bath 
intreaction. In these spin systems typical interaction 
strengths are K\ ~ O.leV and K% < K\, The number of 
nearest neighbours is about 6, depending on the lattice 
structure. Keeping the dominant intreaction (the first 
term) of the qubit with the bath, the problem reduces to 
the effective Hamiltonian given in Eq.2. Consequently, 
our foregoing analysis for the qubit dynamics is applica- 
ble to these spin systems as well. The dynamics in these 
spin systems has been studied at length in the context of 
entanglement generation and propagation[2^. 

We will investigate in Sections III and IV the qubit dy- 
namics exactly for the model Hamiltonian given in Eq.2, 
as a function of the initial nuclear bath spin distribution 



and nuclear polarization strengths in various spin sec- 
tors. As we will see later, through the time evolution, the 
qubit polarization displays a decoherence regime (charac- 
terized by a gaussain decay with a decoherence time scale 
r) for short times. The decoherence time depends on the 
initial bath state. For longer times, the qubit polariza- 
tion revives, as a consequence of the quantum coherent 
evolution of the total system of the qubit and the bath. 
The revival time (Poincare time scale) for the subsystem 
is large in comparison with r. However, for times larger 
than r the sub-dominant interactions, which were ignored 
in the simplified model Hamiltonian, will come into play. 
The revival time becomes larger when the next to leading 
interaction is considered, as the revival time scale is set 
by the smallest interaction strength. Thus the revival 
time in a model Hamiltonian dynamics has no signifi- 
cance for the real system, as the weaker interactions are 
ignored in the model, which will push the revival time fur- 
ther. The decoherence regime boundary, determined by 
the dominant interactions contained in the model Hamil- 
tonian, is not affected too much by the inclusion of the 
sub-dominant interactions. As we will see in Section IV, 
the decoherence time scale from our model Hamiltonian 
compares very well with the experimental numbers, im- 
plying that our model is an excellent approximation. 



III. DECOHERENCE WITH HEISENBERG 
INTERACTIONS 

A. The initial state 

It is straightforward to specify the initial state of the 
central spin and the nuclear spin bath. As mentioned ear- 
lier, we consider the case where a prepolarised electron 
is injected into the quantum dot. Its state is initially 
uncorrelated with the nuclear state. We shall take the 
state of the total system to be a direct product of the 
qubit state and the nuclear bath state. At T = 0, the 
bath will be in its ground state which is sensitive to the 
inter-nuclear interactions. For example, a ferromagnetic 
interction coupling would lead to a fully-polarised bath 
state. This would give rise to a large nuclear magnetic 
field acting on the qubit spin. At nonzero temperatures, 
the bath will be in a mixed state with contributions from 
many bath spin sectors with weights for each spin sector 
Ib (to be fixed by the energetics and the temperature). 
We will consider the case where the intra-bath dynam- 
ics conserves the total bath spin. The density matrix of 
the bath spins may be written, on general considerations 
as pb = ^2^i b PIb where the sum is over all the pos- 
sible values of the total bath spin that the dominantly- 
interacting nuclei can take. In each spin sector the bath 
can have polarizations of various rank. The dynamics 
of the Heisenberg interaction between the qubit and the 
bath, does not mix different bath spin sectors. We can 
determine the dynamics of the qubit and bath in different 
spin sectors separately. 
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FIG. 1: The polarization P A (t) of a spin-1/2 particle A interacting with another spin-1/2 particle B is plotted (left figure) for 
four different initial states of the total system A-B.The initial polarizations of A and B are either 1 or 0.5. The PA(t) and the 
concurrence measure of entanglement, Cab(I), between A and B is plotted (right figure) with time. The initial polarizations 
are given P A (0) = x, P B {0) = (-x + z)/2\/2. 



A further comment is in order here. Consider the nu- 
clei in a state with a given value for the total spin Ib ■ It is 
customary to consider two extreme situations: when the 
nuclear magnetic field has a nonzero expectation value, 
and when only the fluctuating values survive. In the 
latter case the bath is usually dubbed as unpolarized. 
For a maximally-unpolarised bath, the density matrix is 
Pb = jnZ- Then, not only does (Ib) vanish, so do all 
higher order moments involving the multilinears in the 
spin operators. By projection theorem, the expectation 
values of all the moments of the magnetic field operator 
also vanish. There is no a priori nuclear magnetic field, 
and it has to be necessarily induced by the qubit. On the 
other hand, should the system be prepared in such a way 
that (Ib) — 0, but (1^1%) is non vanishing, the system is 
not unpolarized; only the vector polarization is zero. It 
is of course possible that higher rank tensor polarizations 
may survive. We shall clarify the contributions coming 
from tensor polarizations of various ranks, i.e., of various 
higher order magnetic field fluctuations later. 

Finally, we shall take all the nuclei to be spin-1/2 par- 
ticles for simplicity. In reality, the nuclei do have higher 
spins, but the robustness of the analysis is unaffected 
by this approximation. As a warm-up example, let the 
bath be comprised of a single spin- half nucleus, as in the 
following. 



B. Interaction between two spin-half particles 

Though simple, this case is important because it shows 
how entanglement is generated by interactions. The 
Hamiltonian is given by H — KS ■ Ib- Let the initial 



state be a direct product state, 

p(0) = \{X + P A (0) ■ 3 A ) ® Pb(0) -a B ), (5) 

where Pa and Pb denote the initial polarization vectors 
of the qubit and the bath respectively. The polarization 
of the qubit at subsequent times, = Trp(t)aA, is 

easily evaluated to yield, 

P A {t) = cos 2 (Kt/2)P A (0) +sm 2 (Kt/2)P B (0) 

+~sm(Kt)P B (0) x P A (0). (6) 

The polarization for the nuclear spin B is obtained by in- 
terchanging the labels A and B in the above expression. 
Now starting with an initial polarization -Pa(O) = 1, if 
PA(t) < 1 at a later time, the it implies a decoherence. 
However, after decoherence the polarization (of the par- 
ticle A) will build up again, as there is an overall period- 
icity due to the Hamiltonian dynamics. 

It is easy to see from the above equation that PA(t = 
n/K) = Pg(0), and vice versa, which shows that the po- 
larizations get swapped between the qubit and the nu- 
cleus at that time. The decoherence is illustrated in 
Fig. (la) where we have plotted the variation of -Pa(*) 
with time for four different choices of the initial state of 
the total system. Anticipating the analysis in section VI, 
one may note that Pb has a dual role: to cause a pre- 
cession of Pa (represented by the cross product term in 
the above equation) and also change its magnitude. It 
may also be pointed out that since the evolution of the 
total system is unitary, there is a revival of the initial 
state with a period 2tt/K. The "decoherence", however, 
takes place at a much smaller time scale, determined by 
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the first term in Eq|Hl If Pa = Pb, the state remains 
invariant in time. 

Now, the initially-polarized qubit (i.e. in a pure state) 
loses its polarization (i.e. the state is mixed), due to the 
interaction. In other words the total state of the qubit 
and the nuclear spin picks up entanglement. Through the 
time evolution, the entanglement between the two keeps 
changing nonmonotonically. The concurrence measure of 
the entanglement |28[ is easily calculated. In Fig. 1(b) we 
have plotted the entanglement measure, along with the 
polarization Pa(£) as a function of time. From the figure 
it is clear that decoherence sets in (a decrease in the po- 
larization close to zero time) along with the generation of 
entanglement, and at latter times the polarization picks 
up again at the cost of the entanglement. 

The fidelity measure captures how close is pA(t), the 
reduced density matrix of the qubit, to its original state 
Pa(Q)- The average fidelity at time t is given by 

F A {t) = (Tr[ P A(t)pA(0)}/Tr[p 2 A (0)}), (7) 

where the average is indicated over a distribution of ini- 
tial states of the qubit. We have normalized the fidelity 
to unity at time t — 0. If the averaging is done with a 
uniform distribution over all possible initial pure states, 
we have 

F^{t)= l -{l+cos\Kt/2)). 

The average fidelity over all possible initial pure and 
mixed states is obtained as, 

F A {t) = l-i{^- 2 -)s^{Kt/2). (8) 

The decoherence is, of course, not easily apparent in the 
above expressions. The fidelity can be less than unity due 
to either the state is decohered or the state has changed 
coherently (a precession of the polarization vector in an 
external magnetic field). It may not be possible to sepa- 
rate the contributions from a coherent unitary evolution 
and decoherence in the above equation. 

C. Interaction of a qubit with N spin- half nuclei 

The central problem of the paper will be addressed 
here. As already mentioned, we shall take all the nuclei 
to be spin-half, for the sake of simplicity. As argued in 
section II, the Heisenberg interaction involves, not all the 
nuclei, but only those in the immediate neighborhood of 
the localized qubit. 

The Hamiltonian is then simply given by 

H q - h = KS-I B (9) 

The interaction strength K can be written in terms of 
the hyperfine parameters; in general it depends on the 
value of the total spin of the environment. 



Let the initial qubit-bath state be a direct product, 
p(0) = pa{0) <8> Pb(0), where the qubit is denoted by 
A, and the bath by B. The initial state of the qubit 
has the standard form p^(0) = \{X + a a • Pa{0)) where 

Pa(0) = Tr [pa (0)<ta\ ■ The initial bath state ps can be 
written in terms of density matrices corresponding to the 
various bath-spin sectors, ps(0) = Aj B /5/ B (0), which 
displays the initial bath state as an incoherent sum of 
states labelled by bath spin I B , with weights Aj B . The 
nuclear bath state is in general a mixed state in each 
sector of the total bath spin. The bath density matrix 
in a total bath-spin sector Ib can be written in terms of 
tensor polarizations, and spin operators of various ranks 
as given by [30J 

m,n— 1 ' 

(10) 

here, Pj B is the vector polarization of the bath in the 
bath-spin sector Ig, and H" B n is a component of the rank- 
two tensor polarization with the cartesian indices m and 
n. In the above expression we have not shown explicitly 
the terms involving higher-rank tensor polarizations. We 
shall see later that only these polarization terms will be 
relevant to the time evolution of the qubit polarization at 
all times. The components of the second rank tensor op- 
erator are defined by Qf£ = (P B \ l P B :+P 1 ^I B l )/2-I B 5 mn /3 

|30j . The polarizations Pi B and Hi B are determined by 
the expectation values of the vector and tensor spin oper- 
ators respectively in the Is sector, (Ib) — Pi B {Ib + l)/3 
and (Qf£) = n^ n (2I B -l)(2I B + 3)/10. With the above 
definitions the components of both vector and tensor po- 
larization are of order unity. 

The dynamical evolution of the system is governed by 
the equation p(t) = Up(0)U', where the time evolution 
operator is given by (in a sector with the bath spin Ib, 
and setting h — 1) 

U = e- itH <-> =e-^(a lB (t) + b lB (t)S ■ f B )■ (11) 

The time-depedent coefficients are aj B (t) = cos(At) — 
ism(At)/(2I B + 1) and bj B (t) = -Aisia(At)/(2I B + 1), 
where 2A = K(Ib + 1/2). The Heisenberg Hamiltonian 
makes the dynamics exactly solvable. The various com- 
ponents of the total density matrix, corresponding to dif- 
ferent bath spin Ib , can be time evolved separately. After 
determining the state at any time t, a partial trace over 
the bath degrees of freedom yields the expression for the 
reduced density matrix pA(t) of the qubit. We can rep- 
resent the reduced density matrix as 

PA{t) = \{l + PA(t)-S A ). (12) 

The polarization vector Pa (t) carries all the information 
about the qubit, Fa(*) = 1 for a pure state, -Pa(^) < 1 for 
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a mixed state. It depends on the bath-spin distribution, 
and the polarization strengths of the bath in each total 
spin channel. The polarization vector of the qubit at any 
time t is obtained as 



where the gaussian time scale is given by 



P A (t) = f(t)P A (0) + g(t) + h(t) x P A (0) 



(13) 



In the above the sum is over cartesian components 
(m, n = x, y, z) and e n stands for a cartesian unit vector. 
The time-dependent coefficients depend on the bath-spin 
distribution, and are given by 

/(*)= 1 - El B A/. 3(t + l/ + ^ ^ 2 i(i B + |)f ], 
<?(*) = E/ B A /fl 3^±^sin 2 [(/ B + I)f]P jB , 
Ht) = E/ B A/ B 3(f ^+ 1 /2) sin[(/ fl + |)^]P jB , 

n™(*) = E, B A /B H0^sin 2 [(/ B + i)f ]n™. 

(14) 

The coefficient / depends only on the spin distribution, 
and is independent of bath polarizations. All the other 
coefficients depends on various polarizations that may be 
present in the initial nuclear bath state. It is clear from 
EalTSI that P A couples utmost to the second-rank tensor 
polarization. It is a straightforward consequence of the 
Wigner Eckart theorem that the higher rank tensors do 
not couple. We conclude that all nuclear states which 
have the same values of vector and second rank tensor 
polarizations have indistinguishable dynamics. However, 
the dynamics of a central spin-one particle (its state will 
be described by a vector polarization and a rank-2 ten- 
sor polarization) will have contributions from rank-3 and 
rank-4 tensor polarizations of the bath as well. 

For small times the polarization shows a gaussian de- 
cay behaviour (see Fig. 2). This can be seen by expanding 
the time-dependent cofficients shown above. We have 
the leading-order time dependence, / f» 1 — Wft 2 ,g w 
Ky/wht 2 /2, h ps y/w^t, Il mn ps w mn t 2 , where 



w f 



K 2 



= {^Y. X ^{Ib + 1)Pis? 



<W = ^£ \ lB (Ib(Ib + 1) - l)Tl? B n . 



(15) 



(16) 



(17) 



With the above expansion of the time-dependent coeff- 
cients, the small-time behaviour of the polarization shows 
a gaussian decoherence, 



1 _ w h . 2a 

^ = w f ~ sm h 



K^fwh 
2Pa(0) 



pm nil 

E f_A A 
Pa(0) 2 



(19) 



where the angles appearing above are given by, cos Oh = 
h(0) ■ Pa(0),cos6> 9 = y(0) • Pa(0). The decoherence time 
scale r depends on the bath-spin distribution (through 
the moments), and also the magnitudes of the vector and 
tensor polarizations. Nonzero bath polarizations have a 
tendency to increase the time scale over which the gaus- 
sian decay takes place. However, the dominant contribu- 
tion to decoherence comes from the first term (i.e., from 
the function f(t) in Eq.14) which does not depend on 
bath polarizations. 

Another perculiar feature is that the time scale also 
has the initial qubit state dependence through the ap- 
pearance of Pa{0) in the above expression, which is the 
hallmark of nonmarkovian dynamics. Starting with qubit 
pure states (Pa(0) = 1), initially there will be a decay, 
i.e. decoherence, over a time scale r. At later times, 
the polarization will again grow, showing a nonmono- 
tonic behaviour, as we will see in the next section in 
specific examples of the bath-spin distribution. We will 
also consider more general local interations, unlike the 
global interaction considered above, where the qubit can 
interacting with several nuclear spins with the different 
interaction strengths. We will calculate the gaussian de- 
coherence time scale by a direct expansion of the time 
evolution operator in powers of time. 

It should be noted that, for a generic markovian evo- 
lution, an exponential decay is expected. However, the 
case we have considered where the central spin- 1/2 is 
evolved with the spin bath through Hamiltonian dynam- 
ics, a gaussian decay is typical as the process of elimi- 
nation of the bath degrees of freedom is a nonmarkovian 
process. A comparision of markovian and nonmarkovian 
dynamics will be done in the last section, where we con- 
struct a master equation that describes the effective dy- 
namics of the central spin-1/2 particle. 

Before we proceed to discuss special cases and exam- 
ples, we write the expression for the average fidelity for 
the qubit, averaging over all possible initial states of the 
qubit, 



(F A (t)) = l-(~-~)Y XlB f lB[lB 111 sin 2 (J B +l/2)^ 
\ A\ )i v 4 io(J B + 1/2)2 V BT I ) 2 

(20) 



If the averaging is done over all possible initial pure states 
only, we have is confined to the pure states only, 



P A {t) » P A (0)e-^) 2 



(18) 



(F A (t)) = 1 - 1 E ^fsTl^r + 1/2) 



Kt 
~~2 
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IV. EXAMPLES AND SPECIAL CASES 
A. Unpolarized bath 

For an unpolarized bath, polarizations of all ranks are 

identically zero i.e., (I B ) = {Ib f b) = ( J s 7 s-) = °- 
The initial bath state is ps(0) = iprl- For this case, the 
polarisation of the qubit at any time is given by 

P A (t) = f(t)P A (0) (21) 

where f(t) is given in Eq.14 with Xj B — 
^n-Cx T w 7B r +1 ,\ ■ We note that the decay time 

^ "2 IB y+iB + 1 

scale is determined by small time behaviour of f(t). 
Expanding f{t) upto 0(t 2 ) we get 

P A (t) = (1 - NK 2 t 2 /A)P A (0) « e- t2 ^ 2 P A (0) (22) 

where the decay time sclae is r = ■ . . Using the 
known experimental decay time scales of order 10-100 
nanoseconds, and estimating the number of nuclei inter- 
acting with the qubit to be, N ~ 100 (see Section II), 
the typical effective coupling strength relavant for QDQC 
systems works out to be, K ~ 10 _8 eV. These results are 
in broad agreement with the results obtained by |27j . 
This agreement justifies writing an effective Heisenberg 
interaction of the qubit with the total bath spin. 

In Fig. (2) we have plotted P A {t) for an unpolairzed 
bath composed of ]V=10, 30, 50 particles. One can 
clearly see the N dependence on the initial decay. As 
N grows large, the value of P A (t) falls to 1/3 of its initial 
value rapidly, and stays there for a long time. In general 
when Pb(0) — ^\i b Ii b / {21b + 1), the gaussian decay 

time scale is given by Kt = 1/ J < 7| > /3. We can also 

see in Fig. 2 (inset) the revival of the qubit polarization 
over a time scale of the order H/K, which is a consequncc 
of the coherent quantum evolution of the total system 
(the qubit and the bath). However, this revival time 
scale is not a physical time scale for the qubit, as we inl- 
cude the next to leading inteaction strength K' (which 
is exponentially smaller than the dominant interaction 
strength K, see Section II), the revival time will become 
larger ~ h/K'. The decoherence time scale shown above, 
however, will change only slightly due to the inclusion of 
the sub-dominant interactions. Thus though our model 
Hamiltonian is an excellent approximation for the deco- 
herence time scale, the revival times in the dynamics are 
not physical for the real systems. 

In the above examples, the initial bath state is com- 
pletely isotropic, and consequently there is no magnetic 
field produced by the bath. By projection theorem, all 
the moments of the magnetic field also vanish identically 
as argued earlier. We conclude that the decoherence is a 
higher order effect in the following sense: the magnetic 
field produced by the qubit polarizes the bath, which in 
turn produces a magnetic field at the site of the qubit. 
There is no a priori fluctuating magnetic field, contrary 
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FIG. 2: The polarization PA(t) of a spin-1/2 particle, inter- 
acting with a completely unpolarized bath composed of TV = 
10,30,50 spin-1/2 particles, is plotted as a function of time, 
over a complete period. The inset shows the short-time be- 
haviour, a gaussian decay, of the polarization. 

to the statement found in literature. This conclusion 
holds, of course for a larger class of initial nuclear states 
whose vector polarization and rank-2 tensor polarization 
are absent; we have already remarked that the higher 
rank tensors are of no consequence. To illustrate, let 
Pb{0) = J2i B ^i b ^b I + !)• The bath is maximally 
unpolarized in each spin sector. The polarization of the 
bath in a sector Ib can be found easily (by a partial trace 
over the qubit degree of freedom), and the total polar- 
ization of the bath, P B it) = ^ ^Ib Pib > ^ s given by 

(23) 

^From the above expression we can see that the induced 
polarization in each spin sector grows with time; this in 
turn effects the dynamics of the qubit spin. The sub- 
sequent evolution leads to the decoherence in the spin 
state. This is illustrated in Fig. (3), where we have plotted 
the total induced bath polarization, Ps(t), for initially- 
unpolarized nuclear bath with a gaussian bath-spin dis- 
tribution . 



B. Fully-polarized bath 

We now consider the other extreme, a fully polarized 
nuclear bath. All the nuclear spins are parallel. This 
example is of relevance when the inter nuclear coupling 
is ferromagnetic. The resultant magnetic field at the site 
of the qubit can be large, and may be identified with 
the Overhauser field. The semiclassical approach may be 
expected to work the best in this case. We shall present 
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FIG. 3: The nuclear bath vector polarization Ps(f) is plotted 
against time for one complete period, for the case where a 
spin- 1/2 particle is interacting with an initially-unpolarized 
bath composed of iV=100 spin-1/2 particles, and a gaussian 
bath-spin distribution A/ B ~ exp(— 0.1/g 2 ). 

below an exact analysis. The initial bath state here is 
a pure state with bath spin 1b = N/2. ps(0) = | ft 

T)(TT T I - The vector and tensor polarizations 

for this state are 

•W 5/V / _1 °\ 

p N ( ) = _ z;Un(0) = — - 0-10 (24) 

~ y ' N + 2 - V ; 6(iV + 3) I q 2/ 

Substituting these values in Eq^|we get (with A = (N + 
1)K/A), 

IN 

P%{t) = PI(0) + sin 2 (At)(l - 2J*J(0)) 

p . JV_/£+cos2At -sin2Ai \ - 
x W JV + 1 ^ sin 2 At i + cos 2 At J ^ lU j 

(25) 

In the above we have denoted the transverse component 
of Pa by a two-component column vector P±. The ex- 
pressions for the transverse components of the electron 
spin are the same as those obtained by placing the elec- 
tron in a constant magnetic field B — (N + \)K/2z apart 
from terms of order 1/N. For large AT, the polarization 
hardly changes from its initial value. For an inital qubit 
pure state with Pa(0) = —z, we get from Eq.19, the 
gaussian decoherence time scale, r = \[2/ K\[N . These 
results are consistent with the work of Taylor et al[3lj 
(however, there is an apparent discrepancy, which we 
trace to a trivial algebraic error in that paper). 

The reason for such a long-lived polarization of the 
qubit, in the above example, can be traced to the en- 
ergetics of the dynamics. Consider the total state of 
the system to be a direct product of the qubit in down 
spin state and the bath in a fully polarized state. As 



the total z-component of the spin (S z + I z ) is con- 
served in the evolution, the sate at any later time can be 
written as \ip(t)) = ci(t)||)|f ) + c 2 (t)|T)|f - 1), where 
|^r) and j^j— 1) denote the nuclear bath states with 
I z = Yi if — 1 respectively. The energy difference of 
this state with the initial state (with c 2 (0)=0) is given 
by SE/K = N(l - \a{t)\ 2 ) + sfNRe{c* 2 (t)cx{t). Thus, 
for large N, if the coefficient c\(t) changes even slightly, 
there is a large change in the energy, hence the corre- 
sponding transition is not favorable. This implies that 
the qubit polarization cannot change much through the 
evolution. In constrast, if we apply a magnetic field of 
order KN/fiB for the qubit along direction opposite to 
the nuclear polarization, then the qubit polarization is 
given by P\{t) = cos{VNKt)PX(0), which shows that 
the iniial polarization is not preserved. In fact, there will 
be coherent oscillations between the two states shown 
above. This situation can be used in nuclear memory 
where the information coded into the qubit initial state 
at t = can be transfered to the nuclear state through 
the time evolution. These considerations were treated by 
Taylor et alpl]. 
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FIG. 4: The polarization of a spin-1/2 particle interacting 
with a nuclear spin bath with a given total spin, is plotted 
against time. The initial polarization is Pa{0) = —z. The 
two different cases shown are a fully-polarized bath in a pure 
state with the vector and tensor polarizations as shown in 
Eq.25 (solid line), and a bath in a mixed state with a vector 
polarization as shown in Eq.25 and zero tensor polarization 
(dashed line). A nonzero tensor polarization increases the 
initial gaussian decay time scale. 

In the case of a fully polarized bath, the role played by 
the tensor polarization is often glossed over. In fact due 
to the presence of the tensor polarizations (see Eq.25), 
the bath is in a pure state. To highlight the role played 
by the tensor polarization in the above case, and show 
its importance, let us consider the following mixed state 
M0) = ^{1 + ^Pf B I%). For Pf B = ^5 lBtN/2 
this state has a vector polarization same as that of the 
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fully-polarized bath that we considered above, but zero 
tensor polarization. Again considering an initial qubit 
pure state with Pa(0) = —z, from Eq.19 we get the de- 
coherence time scale as r = 2^/3/Ky / N(N + 5) ~ 1/N. 
In Fig. 4 we have plotted Pa(£) when the qubit is interact- 
ing with baths of same vector polarization but different 
tensor polarizations. In the absence of the tensor polar- 
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FIG. 5: The polarization PA(t) of a spin- 1/2 particle inter- 
acting with a bath composed of N = 100 spin-1/2 particles 
and the bath spin distribution \i B ~ exp(— 0.11%). The ini- 
tial polarizations of the bath in various spin sectors are (i) 

Pi B = fii B = o (ii) p Ib = o, n /B / o (iii) p Ib ^ o, fi lB = o 
(iv)p /B + on lB ± o. 



The four cases considered here are (i) no bath polariza- 
tions (ii) no vector polarization but a nonzero tensor po- 
larization (iii) a nonzero vector polarization but no tensor 
polarization (iv) nonzero vector and tensor polarizations. 
The polarizations when nonzero in each bath-spin sector 
are given by 

1 f- 1 °\ 

P, B (0) = z;n lB (0) = - -10 . (26) 

6 \ 2/ 

It can be seen that the decay of the qubit polarization 
is slowest when both vector and tensor polarizations of 
the bath are non-zero, and it is fastest when there are no 
bath polarizations. 
V. SHORT-TIME BEHAVIOUR WITH LOCAL 
INTERACTIONS 

One feature that has emerged from the above is that 
the decoherence (dephasing in the case of ensembles of 
quantum dots) is governed entirely by the short-time 
structure of the evolution operator. This suggests that 
we can use a perturbative solution, to capture the short- 
time behaviour of the polarization. Then we may as well 
enlarge the class of Hamiltonians to include local inter- 
actions and also include an external magnetic field B. 

Consider the Hamiltonian 

JV 

H = S-^Jili+B- S, (27) 



ization one can see that Pa decreases rapidly to zero, and 
in contrast Pa changes very little in the presence of ten- 
sor polarization of the spin bath. The decoherence time 
scale, or the gaussian decay time scale, is much larger 
when tensor polarization is nonzero. 

In Fig. 5, we have shown the short-time behaviour of 
PA(t), with the initial qubit polarization Pa(0) = (x + 
z)/y2 for a few different bath states, in each case with 
a gaussian bath-spin distribution, A/ B ~ exp(— O.L/o). 



which is the most general Hamiltonian describing the hy- 
perfine interactions in quantum dots in the presence of an 
external magnetic field. The exact dynamics of the qubit 
governed by the above Hamiltonian is difficult to solve 
as the time evolution operator does not acquire a simple 
form (given in Eq.ll) as in the case of global isotropic in- 
teraction considered earlier. Taking recourse to a pertur- 
bative approach, we exapnd the time-evolution operator. 
The qubit polarization can be calculated up to 0(t 2 ) as, 



PA(t) = P A (0) - ^{J2MPa(0) x P Bi (0)) + 2P A (0) x B} 

i 

t 2 i ' 

+ j{^2Ji a {Pa i {0)-PA(0)) + J i J j( p A(0) x Pb,{0)) x Pb ( (0) + 2(Pa(0) x B) x B}. (28) 

i i^j 

I 



In writing the above solution we have taken the initial 
state of the system to be a direct product, p(Q) — pa(0)® 
/9b(0), where /9b (0) = pb x ® Pb 2 ' ' ' ® Pb n - The state of 
the i'th nuclear spin is given by pBi — ^{I+&Bi ■-Pb»(0))i 



where its polarization is denoted by Pb ; . 



For an unpolarized bath where all Pb;(0) = 0, the 
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above equation acquires a simpler form 



P A (t) = P A (0)-tP A (0)xg--^2j i 2 P A (Q) (29) 



- ~(P A (0) xB)xB 



^Frorn the above equation one can easily see that P A (t) = 

[1 — x J2i ^ 2 ]-Pa(0), indicating a gaussian decay. More 
interestingly the decay time scale is independent of the 
external magnetic field. 

For the case of a fully-polarized bath where all the 
spins are pointing along the same quantization axis z 
there is no decoherence if the initial qubit polarization 
also points in the same direction, as the state is an eigen 
state of the hamiltonian. On the other hand, if the qubit 
polarization is given by P A {0) = —Z, and the magnetic 
field along any arbitary direction, the polarization can be 
simply read off from Eq. |2HI as 



Pa® = P A (0) - tP A (0) x B 



|^Ji 2 + ^(PA(0)xB)xB 



(30) 



Here the polarization is P A {t) = [1 - ^ Ei ^ 2 ]-Pa(0), 
which is again independent of the external magnetic field. 
From this we can conclude that for short times, the gaus- 
sian decay of the qubit polarization is completely de- 
termined by the sum of squares of the interaction cou- 
pling strengths Jf of the qubit with different nuclear bath 
spins. 



VI. MASTER EQUATION FOR THE 
DECOHERENCE OF THE QUBIT 

In this section we will derive the master equation 
obeyed by the qubit, from the explicit solution for the 
time-dependent density matrix (or the polarization) of 
the qubit given in Eq.13. The advantage in setting up 
a master equation is that it displays the unitary and 
nonunitary parts of the qubit evolution, and their de- 
pendence on the initial bath state, viz. the nuclear spin 
distribution and the polarization strengths. Though the 
master equation carries the same information as the time- 
dependent polarization calcualted above, it displays the 
characteristic effective magnetic field seen by the qubit, 
and the time scale for decoherence and polarization pro- 
cesess explicitly. 

We now proceed to recast our results as the solution 
of a master equation satisfied by the qubit state. The 
most general master e quat ion for a two-level system was 
given by Gorini et.al., |32j and Lindblad [3^. Using the 
notation of Gorini et. al, the master equation for the 



qubit density matrix can be written as 
d 



Of 



PA (t) = -i[H c , PA (t)] 



1 3 

9 ^'{^'mW^i+mm,^]}' 



(31) 



where, the first term represents the unitary component 
of the evolution. The second term is responsible for non- 
unitary processes such as decoherence, polarization and 
equilibration. The coefficients (r^) determine the decay 
or growth of polarization. 

Since we have an explicit solution for the time de- 
pendence of the polarization, the various terms in the 
above equation can be identified. As we will see below, 
the matrix elements Ti m are time dependent, implying a 
non-Markovian evolution, causing a gaussian decay of the 
polarization for small times. Also, the polarization will 
show temporal decay and growth periodically, display- 
ing the underlying Hamiltonian evolution of the spin- 1/2 
and the bath taken together. Let us start by rewriting 
our general solution for the polarization of the spin-1/2 
particle as a function of time given in Ea ll3l as 



P A (t)-g{t)=MP A (0), 



(32) 



where P A , g are the column vectors of P Al g respectively, 
and M is a 3 x 3 matrix given as 



M 



~h z ~\- Tl X y f ~4~ Ilyy h x ~f~ 11^2: 



(33) 



h 



y ~r ±±xz '*>-x ~r ±A -yz f H~ ^z 



IL 



-h-, 



n„ 



Then the equation of motion for the polarization vector 
can be written as 



dP A {t) 
dt 



= DP A (t) + R(t), 



(34) 



where the matrix D = ^fM 1 and the inhomogeneous 
term is given by R(t) — dg(t)/dt — Dg(t). Now, a 
straightforward comparision of this with Eg 13 II yields ex- 
pressions for the matrix elements r; m and the effective 
magnetic field -B e // (in units of h/ns) as 



r; m — Di m + D n 



SlmTrD - i 2J QmnRn, 



and 

mn 



(35) 



(36) 



It is easy to check at small times , that Ti m ~ at and 
\B eff \~[3(l-St2). 

In the above equation for the polarizaion vector, the 
time developement of the different components are cou- 
pled. Let us rewrite the matrix D interms of the sym- 
metric part D s and the antisymmetric part D a as D — 
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D s +D a . The antisymmetric part gives the effective mag- 
netic field, which causes a standard precession of the po- 
larization, with no decay or decoherence. The symmetric 
part causes a decay or growth of the polarization com- 
ponents. We can change to the diagonal basis of the 
matrix D s , and the transformed vectors are P,R. The 
time developement of the polarization components can 
be found directly by integration, and 7,- (eigenvalue of 
D s ) determines the decay or growth of the polarization 
component Pj. For short times, we have — 7$ ~ t, and 
the polarization is given by 



P 2 (t) w ^P?(0)e 2 ^o* 7id *, 



(37) 



which implies a gaussian decay. If all the elements of the 
matrix T (or D) are independent of time, implying time- 
independent decay cosntants 73, then we would have an 
exponential decay which is the hallmark of a markovian 
evolution. 

For an illustration, let us consider a purely vector- 
polarized bath, where the density matrix in each sector of 
the bath spin I B is given by pj B = [X+ j^Pf B I B ]/ + 

1). Here, there is no tensor polarization, and h(t) = h(t)z 
and g(t) — g(t)z. The expressions for the effective mag- 
netic field and the T matrix can be directly calculated for 
this case, 



g fh-fh . 



ifi(f) -I log/ 




(38) 



(39) 



log 



J 



The diagonal elements of the above matrix are related 
to the decay functions through, Tu = 2ji — TrD. We 
have, 




FIG. 6: The decay functions and the effective magnetic field 
are plotted as a function of time when the initial states 
of the qubit and bath are pa(0) = | [)([ | and ps(0) = 
Aj/(2J+ + Pi .z) respectively. The bath is composed 
of N = 100 particles with the bath spin distribution given by 
\i — exp(— 0.17 2 ). Note that, with h = 1, both the decay 
functions 711 and 7x have a dimension of energy. 



of polarization. Though 74 are time dependent, they sat- 
isfy the condition 27j_ > 711 , consistent with the general 
constraint found by Gorini et alj^. 

Finally, we determine the effective magnetic field when 

the initial bath state is given by ps(0) = | ft T)(TT 

t I) an example which we have discussed in section 

IV. B. Here the bath is in a pure state, and the tensor 
polarization is nonzero as discussed in the earlier section. 
In this case, we have 

Beff = {l- ^7 sin 2 (7V + l)Kt/4}™z, (42) 



7i = 72 = 71 = -77 lo g /> 
at 



73-7| l = ^log(/ 2 +/ l 2 ) 



(40) 



(41) 



We see that the polarization of the bath does produce 
an effective magnetic field -B e // which is time depen- 
dent, and causes a precession of PAit). For short times, 
the effective magnetic field is B e ff ps K/^b which is 
about a few gauss using a typical interaction strength 
K ~ 10 _8 eV in QDQC systems. In Fig. 6 we have plot- 
ted the effective magnetic field, and the decay/growth 
functions with time, for a representative gaussian bath- 
spin distribution. As can be seen from the figure, the 
effective magnetic field and 7* show a nonmonotonic be- 
havior, and they take both positive and negative values. 
A negative ji implies a decay of the corresponding polar- 
ization components, and a positive 7; implies a growth 



7|| = 27j 



NK 
(N + l) 



sin(N + l)Kt/2. 



(43) 



Here the strength of the effective magnetic field is about 
a few hundred gauss for N = 100, K ~ 10 _8 eV. 



VII. CONCLUSIONS 

The interaction of a central spin-1/2 particle with the 
nuclear spins in quantum dots is modeled by a Heisen- 
berg exchange interaction. The dynamics of initial direct 
product states of the qubit and the nuclear spin bath 
are investigated. The decoherence of the spin-1/2 par- 
ticle can be seen from the the decay of its polarization 
through the Hamiltonian time evolution. The qubit po- 
larization as function of time is explicitly calculated for 
any nuclear bath with a spin-conserving internal dynam- 
ics, for any bath-spin distribution and any polarizations. 
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The time developement of the qubit polarization has non- 
markovian features. The qubit polarization shows a gaus- 
sian decay /decoherence for short times. For a typical 
interaction strength of the qubit and the nuclear spins, 
K ~ 10 _9 eV, the decoherence time scale is about 100 
nanoseconds. For longer times, the polarization shows 
nonmonotonic behaviour, eventually displaying a peri- 



odicity in time. The gaussian decay time scale depends 
on the bath-spin distribution, a larger width of the dis- 
tribution leading to smaller time scales. The vector and 
tensor spin polarizations that may be present in the nu- 
clear spin bath have a tendency to increase the gaussain 
decay time scale. 
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